Towards a first principles description of phonons in NisoPtso disordered alloys: the 
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Using a combination of density-functional perturbation theory and the itinerant coherent poten- 
tial approximation, we study the effects of atomic relaxation on the inelastic incoherent neutron 
scattering cross sections of disordered NisoPtso alloys. We build on previous work, where empirical 
force constants were adjusted ad hoc to agree with experiment. After first relaxing all structural 
parameters within the local-density approximation for ordered NiPt compounds, density-functional 
perturbation theory is then used to compute phonon spectra, densities of states, and the force con- 
stants. The resulting nearest-neighbor force constants are first compared to those of other ordered 
structures of different stoichiometry, and then used to generate the inelastic scattering cross sections 
within the itinerant coherent potential approximation. We find that structural relaxation substan- 
tially afi'ects the computed force constants and resulting inelastic cross sections, and that the effect 
is much more pronounced in random alloys than in ordered alloys. 

PACS numbers: PACS: 71.20, 71.20c 
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I. INTRODUCTION 

Extensive experimental studies of the lattice dynamics 
of random binary alloys over the past forty yearsiiSiSi^ 
have provided considerable insight into the nature of their 
elementary excitations. Compared with ordered alloys, 
the presence of disorder results in new phenomenona that 
depend on the impurity concentration but also crucially 
on both the relative mass and size difference between 
the constituents. For large mass or size differences, the 
effect of disorder can be dramatic, such as the appearance 
of sharp discontinuities (or split bands) observed in the 
dispersion. 

Among disordered alloys, the NixPti_x system has 
been especially well-studied experimentally^*^ Struc- 
turally simple, the system assumes variants of fee 
over a wide composition range; moreover, the compo- 
nents have large mass (Mpt/MNi=3.323) and size ra- 
tios (rpt/rNi=1.12). Further, the phonons of the end- 
point compounds, elemental metallic Pt and Ni, have 
also been well characterized both experimentally and 
theoreticallyifiiS The interest in these alloys has arisen 
in part because homogeneous crystals are easy to grow at 
nearly any concentration, making the system attractive 
for experimental study. 

Theoretically, there have been many attempts to de- 
scribe the complex nature of phonon excitations in ran- 
dom alloys. The majority were based on the coher- 
ent potential approximation (CPA)2i and its various 
generalizationspi^Jii^iiiiiiSiiSiii The CPA is a single- 
site, mean-field theory capable of dealing only with 
mass disorder; it is commonly generalized to treat both 
ofT-diagonal and environmental disorder. A transla- 
tionally invariant CPA based on the augmented space 
formalism^^ has recently been used successfully to de- 
scribe the lattice dynamics of Ni55Pd45 and NisoPtso 
random alloysji^ This formalism, known as the itiner- 



ant CPA (or ICPA), captures the effects of both force 
constant and mass disorder. 

A key deficiency of the application of those theories is 
that the Green's functions are constructed from empiri- 
cal force-constants. Moreover, in Ref. |l9| empirical force 
constants were then themselves adjusted to fit frequen- 
cies and linewidths extracted from neutron-scattering 
data. A better alternative to reliance on phcnomcnolog- 
ical force constants is to compute them from first prin- 
ciples. In recent years, first-principles density-functional 
theory (DFT) has achieved the ability to predict accu- 
rately material properties without experimental input 
or adjustable parameters. A wealth of first-principles 
studies has demonstrated that electronic, vibrational, 
and transport properties are often extremely sensitive to 
the precise details of the atomic arrangement iSSiSiiSSiSiS^ 
Providing a parameter-free description of the vibrational 
properties of binary alloys would allow a deeper, atomic- 
scale understanding of the phenomenology. 

An especially relevant atomic-scale aspect of these dis- 
ordered alloys which can be captured by first-principles 
methods is the degree to which the atoms displace from 
their high-symmetry fee sites because of size mismatch; 
these atomic "relaxations" would be expected to have 
considerable impact on the force constants. Experimen- 
tally, non-negligible displacements are routinely found in 
binary alloys where the size difference between the con- 
stituents is large. For example, EXAFS measurements 
of CuPd alloys have shown that 1 at% Pd in Cu changes 
the Cu-Pd nearest neighbor distance to 2.560 A, from 
the 2.515 A Cu-Cu bond distance in pure Cu^ Another 
EXAFS study of the Ni^Aui-x system revealed three dis- 
tinct Ni-Ni, Au-Au and Ni-Au bond lengthsiS& 

Within the first-principles framework, some degree of 
translational symmetry is almost always assumed and 
disordered systems, such as random binary alloys, are 
treated approximately, usually by averaging properties 
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over several different ordered atomic configurations. Cur- 
rent computational capabilities limit the size of each con- 
figuration, or supercell, and also the number of configu- 
rations sampled. However, information from one partic- 
ular ordered compound can nonetheless provide insight 
into the physics of its disordered counterparts. Although 
the force constants of a particular ordered system are 
not necessarily expected to be transferrable to a random 
environment, they should yield insight into the relative 
contribution of lattice composition and relaxation to the 
vibrational properties of NisoPtso, a fundamental feature 
of these systems which has yet to be studied. 

In this article, we attempt to capture the effects of 
atomic relaxation on the inelastic incoherent neutron 
scattering cross section of disordered NisoPtso using first- 
principles force constants computed from ordered alloys 
as input to the ICPA. We build on previous work with the 
ICPA in which empirical force constants were adjusted in 
an ad-hoc manner to agree with experimental cross sec- 
tion spectrapiS With a large mass and size ratio, NisoPtso 
is expected to possess substantial atomic relaxation and 
large force constant ratios; further, we have a wealth of 
experimental and theoretical results for comparison. 

The paper is organized as follows. In section II, we 
briefly describe the theoretical tools used in this work to 
connect the first-principles calculations with the random 
alloy systems. In Section III, we present phonon spec- 
tra, site-projected phonon densities of states, and force 
constant results for the ordered NisgPtso compounds in 
pseudo-cubic and LIO structures; the force constants of 
Ni25Pt75 and Ni75Pt25 compounds in L12 structures are 
calculated for comparison. We then investigate the effi- 
cacy of ordered-alloy force constants in reproducing the 
spectrum of the random 50-50 alloy. Concluding remarks 
appear in Section IV. 



II. THEORETICAL BACKGROUND 

A. Details of first-principles calculations 

In what follows we report computations of the ground- 
state properties and phonon spectra of NisoPtso ordered 
alloys in the LIO structure. We consider two cases. In the 
first, the c/a ratio of the tetragonal LIO structure is set 
equal to unity and only the volume is relaxed. We have 
referred to this structure above as "pseudo-cubic", and 
below it is simply referred to as "cubic" . In the second, 
the c/a ratio is relaxed as well. Structural information 
and force constant data for pure Ni and Pt and ordered 
NiaPt and NiPta are also presented. We use density- 
functional theory within the local-density approximation 
(LDA) to relax the different structures. We employ a 
plane-wave pseudopotential approach with the Perdew- 
Zunger parametrization^ of the LDA as implemented in 
the PWSCF packagei2£ Ultrasoft pseudopotcntials^^ arc 
used for Ni and Pt and explicitly treat 10 valence elec- 
trons for each species with a kinetic energy cutoff of 30 



Ry. Nonlinear core corrections are included in the Ni 
pseudopotentiali21i The Brillouin zone (BZ) integrations 
are carried out with Methfessel-Paxton smearing^! us- 
ing an 8x8x8 k-point mesh, which corresponds to 70 k- 
points in the irreducible wedge. The value of the smear- 
ing parameter is 0.03 Ry. Hellmann-Feynman forces are 
calculated, and the atoms arc relaxed steadily toward 
their equilibrium values until the forces are less than 1 
mRy/au. These parameters are found to yield phonon 
frequencies converged to within 5 cm~^. 

Once adequate convergence is achieved for the ground- 
state structural properties, the phonon spectra are ob- 
tained from linear response using density-functional per- 
turbation theory (DFPT)i^ Within DFPT, the force 
constants are conveniently computed in reciprocal space 
on a finite q-point grid; to obtain the real-space force- 
constants, Fourier transforms are then performed numer- 
ically. The number of unique real-space force constants 
and their reliability depend on the density of the q-point 
grid: the closer the q-points are spaced, the more ac- 
curate the force-constants of the higher neighbors will 
be. In this work, the dynamical matrix is computed on 
a 6x6x6 q-point mesh commensurate with the k-point 
mesh for all structures i2i 



B. Itinerant Coherent Potential Approximation 

The itinerary coherent potential approximation 
(ICPA) is a Green's-function based technique for treat- 
ing random substitutional alloys. It is capable of treat- 
ing off-diagonal and environmental disorder and there- 
fore appropriate for the study of phonons in alloys with 
strong force-constant disorder. The ICPA was developed 
and used to study NWPd45 and NisoPtso alloys. The 
recent results of Ref. [ly demonstrate its advantages over 
the simpler CPA for understanding both the dispersion 
and lifetimes of phonons in random binary alloys. How- 
ever, in that work empirical force constants were adjusted 
to fit frequencies and linewidths extracted from neutron- 
scattering data. Using DFPT to calculate the phonon 
frequencies and force constants of various ordered struc- 
tures allows us to make a series of comparisons with both 
the empirical force-constant results and with experiment. 
In doing so, we test the extent to which ordered-alloy 
force constants resemble random-alloy force constants 
and thereby illuminate the effect of relaxation on force 
constants, as reported in the next section. 



III. RESULTS AND DISCUSSIONS 

A. Structural information and ground state 
properties 

Solid elemental nickel and platinum stabilize in the 
fee structure under standard conditions and are mutu- 
ally soluble at all concentrations. When combined, they 
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FIG. 1: Schematic L12 (left) and "cubic" (right) structures. 
For NiaPt the white circles in the L12 structure indicate Pt 
atoms, the gray ones Ni, and vice versa for NiPta. The white 
circles in the "cubic" structure indicate Pt atoms, the gray 
ones Ni atoms. In the L12 structure, atoms of each species 
form three dimensional next-neighbor networks; the "cubic", 
on the other hand, is a layered structure of Ni-Pt layers in the 
(001) direction with c = a. In the LIO structure c/a is free to 
relax, allowing an adjustment of the Ni-Pt layer separation. 



form ordered compounds at 1:3, 1:1 and 3:1 Ni:Pt atomic 
proportions. The 1:3 and 3:1 phases crystaUize in the 
cubic L12 (CusAu) structure; the 1:1 phase takes up the 
LIO crystal structure24 as shown in Fig.l. All of these 
structures contain 4 atoms per unit cell. Their electronic 
structures, ordering tendencies, and magnetic properties 
have been investigated in detail over the years. '^^ 

In Table I, we report the equilibrium lattice parame- 
ters for the three compounds obtained from first princi- 
ples calculations performed as described above and com- 
pare them with the experimental datai^ The LDA lattice 
parameters are 1 to 2 % less than experiment, which is 
within the normal range of error of the LDA and qualifies 
as very satisfactory agreement. 



TABLE I: Lattice parameters of NiPt ordered compounds cal- 
culated using DFPT and compared to experimental values 
taken from Ref. ,36- "Cubic" refers to the lattice constant a 
of the NiPt structure obtained by fixing c/a at unity. 



System 


a (au' 


Expt. c/a Expt. 


Ni 


6.480 


6.650 


NiaPt 


6.744 


6.892 


LIO NiPt 


7.136 


7.244 0.934 0.939 


"cubic" NiPt 


6.99 




NiPta 


7.200 


7.320 


Pt 


7.394 


7.410 




B. 


Phonon spectra 



The phonon spectrum of the 50-50 ordered alloy in the 
LIO structure consists of twelve branches, three acoustic 
and nine optical. The results of our first-principles DFPT 
calculations are shown along the high symmetry lines of 
the simple cubic BZ in the top panel of Fig. 2 for the 
"cubic" structure, i.e. LIO with c = a. Results for the 



fully relaxed LIO structure are shown in the bottom panel 
of Fig. 2. The corresponding densities of states (total and 
component-projected) are shown on the right hand side 
of each figure. 

Although many features of the "cubic" and LIO dis- 
persion curves are similar, several key differences emerge 
upon close examination. The frequencies of the optical 
branches in all three principal directions F-X, F-M and F- 
R are higher in LIO than in the "cubic" arrangement. In- 
deed, since the unrelaxed "cubic" structure has a slightly 
larger volume than the LID (by roughly 0.5%), it would 
be expected to possesses frequencies that are, in an aver- 
age sense, smaller than those seen in the LIO. Comparing 
the eigenvectors and the projected densities of states of 
the two structures, we find that, for both structures, the 
longitudinal optical (LO) branches are heavily dominated 
by the vibrations of the lighter Ni atoms, and the acoustic 
modes are rich in the motions of heavier Pt atoms, again 
as expected. However, closer scrutiny of the LIO dis- 
persion curves reveals a much more appreciable Pt con- 
tribution at intermediate frequencies ('~ 4-5 THz). The 
enhanced Ni-Pt mixing at these frequencies is consistent 
with the fact that in LIO, the interplanar distance de- 
creases from 4.9427 to 3.335 au. 

In the "cubic" structure, component-projected densi- 
ties of states indicate low frequency phonons are domi- 
nated by Pt and the high frequency phonons are domi- 
nated by Ni, with very little evidence of available states 
around the middle. This stricter separation of the Ni 
and Pt contributions is the signature of split-band be- 
havior. Since no such split-band behavior is seen in ran- 
dom NisoPtso alloyr^ the force constants computed from 
the "cubic" structure would not be expected to represent 
those of the random alloy. Relaxing to LIO, on the other 
hand, results in contributions from both Ni and Pt in 
the low-frequency acoustic branches across the Brilloun 
zone, consistent with the appearance of more states in the 
middle of the spectrum; the branches are more evenly dis- 
tributed in energy. Similarly, the optical branches, also 
very closely spaced in the "cubic" structure, are seen to 
broaden in LIO. These results illustrate how structural 
relaxation modifies the vibrational spectra of the ordered 
alloy, bringing it closer to that of the disordered alloy 
within which atomic relaxation is less constrained. Table 
II lists the computed real-space force constants for the 
elements, the ordered phases, and, for comparison, the 
empirical force constants of Ref. for the disordered SO- 
SO alloy. We first consider the composition dependence 
of the force constants of the fully relaxed structures in 
the sequence Pt, NiPtg, LIO NiPt, NigPt and Ni. All 
Pt-Pt and Ni-Ni force constants increase monotonically 
in magnitude with Ni concentration, the Pt-Pt and Ni- 
Ni nearest neighbor separations decrease concomitantly, 
compressing the atoms, roughly speaking, and causing 
the force-constant increase. The Pt-Pt force constants, 
however, increase substantially less than the Ni-Ni force 
constants. The ratio of maximum to minimum force con- 
stants ranges from 1.24 (zz-component) to 1.47 (xx) for 



FIG. 2: Dispersion curves and phonon densities of states for "cubic" NiPt (top) and relaxed NiPt in the LIO structure 
(bottom), calculated from DFPT within the LDA. The leftmost panels are the dispersion curves along high-symmetry points 
in the Brillouin zone. The panels on the right are the total and component-projected densities of states. In the panels showing 
the dispersion curves, the solid lines are the transverse acoustic modes, the dotted lines are the longitudinal acoustic modes, 
the dot-dashed lines are the longitudinal optical modes, and the long-dashed lines are the transverse optical modes. In the 
panels showing the densities of states, the dashed lines are the Ni contributions, the dotted lines the Pt contributions, and the 
solid lines are the total contributions to the densities of states. 



TABLE II: Real-space nearest neighbor force constants for the elements and the three compounds calculated with DFPT; the 
empirical (emp) force constants of Ref. for the random NiPt alloy are included for comparison. The units are dyn cm~^. 





Pt 


NiPta 


"cubic" NiPt 


LIO NiPt 


NiPt (emp) 


NiaPt 


Ni 




Ni-Ni 






-4863 


-5252 


-15587 


-11758 


-13800 


Ixx 


Ni-Pt 




-13200 


-17418 


-13279 


-13855 


-29762 




Ixx 


Pt-Pt 


-27744 


-37571 


-51175 


-40158 


-28993 






Ixx 


Ni-Ni 






318 


498 


436 


253 


149 


Izz 


Ni-Pt 




2813 


2171 


2643 


348 


3044 




Izz 


Pt-Pt 


5512 


5559 


9486 


6623 


7040 






Izz 


Ni-Ni 






-5373 


-6801 


-19100 


-13028 


-15530 


Ixy 


Ni-Pt 




-15949 


-21885 


-15859 


-15280 


-35310 




Ixy 


Pt-Pt 


-30969 


-40474 


-59322 


-45529 


-30317 






Ixy 



Pt-Pt, whereas that for Ni-Ni ranges from 2.28 {xy) to 
3.34 {zz). We note also that the composition dependence 
of the Ni-Pt force constants is weakly nonmonotonic for 
the zz and xy components with little difference in all 



three force constants between their values from NiPta 
and LIO NiPt. 

In Ref. 0, an excellent fit to the experimental coher- 
ent and incoherent inelastic scattering crossections at 
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the 50-50 composition was obtained through adjustment 
of the measured Ni-Ni, Ni-Pt, and Pt-Pt nn force 
constants. In particular, the Ni-Ni nn force constants 
had to be reduced well below their experimental values 
in pure Ni to obtain good agreement, while adjusting the 
Pt-Pt nn force constants was found to be substantially 
less effective in modifying the cross sections. The Pt-Pt 
force constants were thus modified only slightly from 
their experimental values in pure Pt, as seen in Table 
II. (Details of the fitting strategy are given in 
This receives some post hoc justification given the 
relative insensitivity of the Pt-Pt nn force constants to 
composition manifest in Table II and noted above. We 
comment further in the next section on the implications 
of the sequence of force constants for the three 50-50 
compositions in the table. 



C. Ordered alloy force constants as a first 
approximation to the disordered alloy 

A detailed understanding of phonon excitations in dis- 
ordered alloys has been impeded by absence of detailed 
knowledge of the force constants. Even in the case of 
substitutional disorder, the atoms have the freedom to 
relax, and, due to the random occupation of sites, the in- 
teractions among various species are expected to be quite 
different than in the ordered alloy. 

To explore the impact of lattice relaxation on the 
force constants, we calculate the incoherent inelastic 
scattering crossection of the random NiPt alloy using 
the first-principles force constants of the unrelaxed "cu- 
bic" structure and the fully relaxed LIO structure. The 
results are shown in Fig. 3. Note that the "cubic" curve 
is considerably higher than the data at low frequencies, 
and lower at high frequencies. These discrepancies are 
substantially reduced by relaxation: the force constants 
from LIO markedly improve the agreement between ex- 
periment and theory. Detailed differences do remain be- 
tween the LIO and the experimental curves, however. 
In particular, the cross section is still higher at than 
the data at low frequencies (below about 4 THz), and 
the high frequency cutoff predicted from the LIO force 
constants is lower than measured, a disagreement un- 
changed by relaxation. We attribute these discrepancies 
to the constraints on relaxation enforced by the order 
on the nearest neighbor separations. In the random alloy 
such constraints are obviously absent, and force constants 
computed from larger supercells averaged over several 
configurations may improve the agreement. For exam- 
ple, allowing the Ni-Ni nn separation to drop below that 
predicted from our ordered LIO structure would increase 
the Ni-Ni force constants, shifting spectral weight from 
lower to higher frequency in Fig. 3 and improving the 
agreement with the empirical force constants in Table II. 
Similarly, the Pt-Pt nn separation can increase, which 
would tend to reduce the Pt-Pt force constants toward 



their empirical values in Table III. The case of the Ni-Pt 
force constants is more subtle. The Pt-Pt force constants 
are larger in magnitude than the Ni-Ni force constants, 
which suggests that the energetic cost of compressing Pt 
is larger than that of expanding Ni. This suggests in turn 
that the Ni-Pt separation increases in the random alloy, 
which would also reduce the Pt-Ni force constants, as we 
find. 

Several other interesting trends can be seen in the force 
constants. The lattice constant of the "cubic" NiPt alloy 
is 6.99 a.u. Due to the contraction of the "cubic" lattice 
relative to the tetragonal LIO structure, the larger Pt 
atoms have less space to vibrate, resulting in a stiffening 
of the Pt-Pt interaction, as can be seen in Table II. In 
contrast, since Ni is a relatively much smaller ion than 
Pt, the Ni-Ni interaction does not significantly change. 
For this reason we predict too much weight in the low- 
frequency region of the scattering cross section in Fig. 3. 
The low-frequency region is dominated by Pt atoms and, 
due to the considerable hardening of Pt-Pt compared to 
Ni-Ni, the Pt-Pt interaction dominates throughout the 
entire frequency spectrum for "cubic" NiPt. After relax- 
ation to the LIO structure, the situation improves. The 
atomic separation in the Ni and Pt planes grows (since 
the in-plane lattice parameter increases), softening the 
Pt-Pt interactions and producing better agreement with 
the experiment. These results underline the importance 
of capturing relaxation effects in random environments. 

While we have captured some of the relaxation occur- 
ing in the random alloy, a deficiency of our theory is cer- 
tainly that we fail to capture all of it. For example, we 
neglect randomness in the relaxation , and this may be 
responsible for our underestimate of the position of the 
high-frequency band-edge in Fig. 3. As mentioned, the 
Ni-Ni interaction is still too weak to push the upper band- 
edge towards higher frequency. It also has to be noted 
that the shapes and the features of the two scattering 
cross section curves obtained by ICPA using the "cubic" 
and LIO force constants have very little similarity, imply- 
ing that structural relaxation is even more pronounced in 
random alloys than in ordered systems as stated above. 
The incoherent scattering cross sections are nothing but 
the component-projected densities of states weighted by 
their incoherent scattering lengths and therefore carry 
the signature of the component-projected densities of 
states. Fig. 2 suggests that though quite dissimilar in 
their various features, there is more overall similarity in 
the densities of states curves for the ordered structures 
than to that of their random counterpart. 

Since the use of LIO ordered-alloy force constants im- 
proves agreement with experiment, it is tempting to fol- 
low the same procedure for NixPti_x at a;=0.25 and 
a;=0.70, other concentrations at which experiments have 
been done. However, in the simplest ordered counterpart 
of Ni25Pt75, a Ni atom does not have any Ni nearest- 
neighbors and the Ni-Ni force constant has to be esti- 
mated empirically. Moreover, Ni7oPt3o does not even 
have a simple ordered counterpart. Thus the agreement 




FIG. 3: Inelastic incoherent scattering cross sections for NisoPtso random alloy. The solid line is the experimental result, the 
bold dashed line is the ICPA result with empirical force constants, the dashed line is the ICPA result with first-principles force 
constants for LIO NiPt and the dotted line is the ICPA results with first principles force constants for "cubic" NiPt. 



for the 50-50 case may be of limited general use. The 
force constant data in Table II, however, suggest that 
the Ni-Ni and Pt-Pt force constants vary approximately 
linearly with composition, following a Vegard's law of 
sorts for force constants. Assuming this holds through- 
out the concentration range, one could then easily in- 
terpolate to an arbitrary concentration and use the re- 
sults as a reasonable first guess for force constants at that 
particular concentration. On the other hand, the Ni-Pt 
force constants are approximately the same for NiPta and 
NiPt(LlO), justifying a linear interpolation for starting 
values for NiaoPtyo alloy. 

As alluded to above, another route would be to con- 
sider more complicated ordered structures. Zunger and 
collaborators have demonstrated that a clever choice of 



supercell, plus additional relaxation can reliably repro- 
duce optical and thermodynamic properties of certain 
binary alloys n^^i^Z, Their supercell was constructed by 
arranging a minimal set of atoms so that the first few 
peaks of the radial correlation function matches that of 
the disordered alloy for a given concentration. Despite 
a plethora of works successfully using these and simpler 
atomic arrangements to model the properties of alloys, 
to our knowledge there has yet to be a study assessing 
the ability of any of these structures to reproduce the 
details of the phonon spectra of random alloys. In a sub- 
sequent communication, we shall examine the utility of 
using such an artifical supercell to obtain force constants 
in a random environment. 



IV. CONCLUSIONS 



In a first attempt to reach the greater goal of obtain- 
ing reliable first- principles force constants for NisoPtso 
alloys, we have studied the vibrational properties of or- 
dered stoichiometric compounds of Ni and Pt using first 
principles DFPT. The analysis of the phonon spectra 
based on the interactions between different constituents 
of the alloy and their variations under different environ- 
ments provides useful insight about those interactions in 
random environments. The calculated force constants 
of compounds of different stoichiometry elucidate two 
other important features of the force constants: first, 
they provide a posteriori validation of our choice of phe- 
nomenological force constantglS, for NispPtso; and sec- 
ond, they clarify the composition-dependent variation 
of the force constants of ordered, stochiometric com- 
pounds. This paves the way to reasonable estimates 
of off-stoichiometric force constants at arbitrary concen- 
tration. Most importantly, this study strongly suggests 
that a successful theoretical description of the vibrational 
properties of random alloys must include the effects of 
atomic relaxation. 
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